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SIGNIFICANCE  AND  EXPLANATION 


•w 


The  following  result  is  useful  in  connection  with  matrix  applications:  If 
A^ (A)  denotes  the  largest  eigenvalue  of  a hermitian  matrix  A,  then 

(1)  A^A+B)  £ A1(A)  + A^  (B)  , 

i.e.,  if  A and  B are  hermitian  matrices,  the  largest  eigenvalue  of  A+B  is 
at  most  the  sum  of  the  largest  eigenvalue  of  A and  B.  The  quantity  A^(A)  is 
a functional,  i.e.  a scalar  depending  on  the  matrix  A.  The  above  example  suggests 
the  following  problem  which  is  solved  in  this  paper:  Determine  all  functionals 
41(A)  depending  only  on  the  eigenvalues  A^,...,!^  of  A such  that  4>(A)  is 
convex , i.e. 

(2)  4>(aA  + (l-a)B)  £ a 4>(A)  + (l-a)4>(B)  , 0 £ a £ 1 

when  A,  B are  hermitian.  In  economics  and  biology  one  very  often  deals  with 
non-negative  matrices.  Denote  by  A^fA)  the  spectral  radius  of  a non-negative 
matrix  A £ 0,  i.e.  the  largest  non-negative  eigenvalue  of  A.  The  fact  that 
A^A)  > 1 or  A^(A)  < 1 plays  a crucial  role  in  the  stability  behaviour  of  the 
system.  So  any  convexity  results  on  A^(A)  are  helpful  to  estimate  A^(A).  Un- 
fortunately (1)  does  not  hold  in  general  for  A,  B non- negative.  In  this  paper 
we  prove  the  validity  of  (1)  for  A,  B non-negative  if  B-A  is  a diagonal  matrix. 
We  extend  this  result  for  more  special  type  of  non-negative  matrices.  To  derive 
these  results  we  bring  new  characterizations  of  the  spectral  radius  of  non-negative 
matrices. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC,  and  not  with  the  author  of  this  report. 
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the  following  order 


Thus  X (A)  is  a convex  functional  on  H . Ky  Fan  extended  (1.2)  [3] 


max  l 

(x  ,x.)=6..  i=l 


is  a convex  functional  on  H . A function 


is  called  a spectral  function  if 


is  a convex  functional  on  A.  It  turns  out  that  F must 


we  characterize  all  F for  which 
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(1.8) 


n n 

l «4  ■ l h 

i=l  i-1 


We  also  characterize  all  F such  that  $ is  strictly  convex.  Let  A be  an  n x n non- 
negative matrix.  As  usual  denote  by  r(A)  the  spectral  radius  of  A.  So  X^A)  - r(A)  . r(A) 
is  not  a convex  functional  on  non-negative  matrices.  For  example  consider 

U-9)  A = { “ t ) , r (A)  = /e  . 


:)• 


Recently  [1]  Cohen  proved  that  r(A)  is  a convex  function  in  i-th  diagonal  entry  of  A for 

any  1 < i £n.  We  extend  Cohen’s  result  namely,  we  show  that  r(A+D)  is  convex  on  D — 

n 

the  set  of  all  n x n real  diagonal  matrices.  In  fact  this  result  is  a conseouence  of  the 
Donsker-Varadhan  characterization  of  r(A)  [2] , In  Section  3 we  bring  more  general  characteri- 
zations of  r(A)  by  using  a certain  fundamental  inequality  for  non-negative  matrices  establish- 
ed in  [5) . This  enables  us  to  show  that  log  r(eDA)  is  also  convex  on  D for  a non-negative 
A.  If  A happens  to  be  an  M-matrix  then  we  have  a stronger  result.  Namely,  r (DA)  is  con- 
vex on  D+  — the  subset  of  non-negative  matrices  in  D . This  is  done  in  Section  4. 
n n 

In  Section  5 we  show  how  the  results  of  Section  2 can  be  extended  to  the  non-symmetric  case 
by  assuming  that  A is  a totally  positive  matrix  of  order  j(TP^).  We  shall  state  our  results 
in  case  that  A is  a TP  (=TP  ) matrix.  That  is  all  minors  of  A (of  all  orders)  are  non- 
negative. In  that  case  we  have 

(1.10) 

If  A is  non-singular  then  the  last  ineouality  is  strict.  Let 

(1-11)  *<D)  = F ( log  X (e°A) log  X (eDA>)  . 

l n 


X,(eDA)  > X- (e°A)  >...>  X (eDA)  > 0,  D e D 
i — z — — n — n 


Then  <J>  is  convex  on  A c if  and  only  if  F is  convex  on  X and  Schur’s  order  preserving. 

We  remark  that  the  results  in  Section  2 hold  for  symmetric  compact  operators  in  Hilbert 
space.  The  results  of  Section  3-5  can  be  extended  to  appropriate  integral  operators,  for 
example,  as  it  was  pointed  out  in  [5] . 
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2 . Con-'ex  functions  on  the  spectrum  of  hermitian  matrices 


Let  A be  an  n » n hermitian  matrix.  We  can  view  A as  a self  adjoint  operator  on  <t 


endowed  with  the  standard  inner  product 


# n 

(x,y)  = y x , x,y  e a 


Since  the  eigenvalues  of  A sure  real  we  arrange  them  in  the  decreasing  order 

(2.2)  X. (A)  >...>  X (A)  . 

l — — n 

Denote  by  £, the  corresponding  set  of  orthonormal  eigen-vectors  of  A 
l n 

(2.3)  A£a  = Xi(A)Ci,  = «i;.,  i,j  = l.....n  . 

Let  H denote  the  set  of  all  n x n hermitian  matrices.  Since  X, (A)  has  the  maximal 
n 1 

characterization 

X (A)  = max  (Ax,x)  , 

(x,x)=l 


X, (A)  is  a convex  function  on  H . More  generally  we  have  [4] 

1 n 

Theorem  2.1.  Let  {c^}"  be  a decreasing  sequence  of  real  numbers 


a,  > a.  >...>  a . 
1 — 2 — — n 


Then  for  any  A belonging  to  H 


1 a X (A) 
i=l 


(x.,xj)-6.ji,j-l,...,n  i 1 


£ ai(Axi,xi> 


Assume  that  the  equality  sign  holds  for  some  x^,...,x  . Let 

(2.6)  «!-•••-  > «ii+1  ="•'  “i2  > — > \ * Vli0=0)  • 

Then  there  exists  an  orthonormal  eigensvstem  of  A such  that  the  following  subspaces  coincide 


(2.7)  U ) = [x  ,...,x  1,  j = 0,...,r-l  . 

V Xj+1  Xj  1 3+1 

The  characterization  (2.7)  in  the  case  that  =...=  ou  = 1,  <»i+1  =...=  an  = 0 was  established 


by  Fan  (31  . 


“ * iwypti  ■ 


In  particular 


(2.8) 


n 

♦(A)  = l a A (A) 
i-1 


is  a convex  functional  on  H if  (2.4)  is  satisfied.  That  is 

n 


(2.9)  ♦ (cA  + (l-c)B)  < c*(A)  + (l-c)<J>(B),  A,  B H , 0 < c < 1 . 

— n — — 

We  now  are  ready  to  state  the  problem  which  we  solve  in  this  section.  A function 

(2.10)  * : A A <z  H 

n 


is  called  a spectral  function  if 

(2.11)  4(A)  - F( X (A) ,...,A  (A))  . 

l n 


That  is  <t>  is  defined  on  the  spectrum  of  A.  Our  problem  is  to  characterize  all  convex 
spectral  functions  on  H^.  To  answer  this  problem  we  introduce  some  notation  and  definitions 
Let  a = and  6 - (8j*...#8n)  be  two  vectors  satisfying  (2.4).  According  to 

[7,  sec.  2.18]  a is  majorized  by  6/  which  is  denoted  by  a < B,  if 

k k 

(2.12)  J a.  <_  J 8.,  k = l,...,n-l  , 

i-1  i-1  1 


(2.13) 


Denote 

(2.14) 


n n 

I l 

i-1  i=l 


A (A)  = (X  (A)  , . . . ,A  (A) ) 
1 n 


From  Theorem  2.1  we  obtain 

Lemma  2.1.  Let  A , B £ H . Then 
n 

(2.15)  X (A+B)  < X (A)  + X (B) 
Moreover . 

(2.16)  X ( A+B ) = X (A)  + X (B) 
if  and  only  A and  B have  a common  eigenvector  system 


| - -.v 
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t' 


(2.17) 


A?i  - Xi(A)ei,  B?i  - Xi(B)Cl,  U^)  - «jy  i.j  - lf*«n  . 


t 


Proof.  Let 


(2.18) 


(A+B)  ■ X^(A+B)(^i  (C^/£^)  " j i i*j  — 1»  • • • #n 

So  for  any  a = (a^,...,<0  which  satisfies  (2.4)  we  get 


(2.19) 


n n 

l a X (A+B)  - l a ((A+B) 5., e.)  1 
i-1  11  i-1 

n n 

l a X (A)  + l a X (B)  . 

i-1  1 i-1  1 1 


This  establishes  (2.15).  Suppose  that  (2.16)  holds.  Then  we  must  have 
n 


(2.20) 


I °iXi(A)  “ S ai(A5i'5i)'  I aiXi(B)  - I cti<Bei»Ci) 

i-1  l-l  i=l  i=l 


for  any  a 2 ^...>  a^.  Choose  = n-i.  Then  the  equalities  (2.7)  imply  (2.17).  This 

conclusion  is  in  fact  is  stated  in  Theorem  3.1  in  [4]. 

By  k"  denote  the  following  subset  of  »n 


(2.21) 

Clearly 

(2.22) 

Let 

(2.23) 


*>  ={*!*=  <*i V'  X1  ^X2  i"'^Xn} 


X : H -*•  K"  ( X (A)  = (X,(A),...,X  (A)))  . 

n > x n 


\U)  = x . 


Thus  the  function  F in  terms  of  which  4>  is  constructed  satisfies  F : X -*•  R . 

Let  D be  the  set  of  all  n x n real  diagonal  matrices  and  D X the  set  of  all  diagonal 


matrices 

(2.24) 


D(a)  = diag  {a,,...,a  },  a,  > a„  >...>  a 
l n 1 — 2 — — n 


Suppose  0 £ X.  Then  D(B)  e A.  Thus  the  assumption  that  4>  is  convex  on  A implies  in 

particular  that  4>  is  convex  on  D*  n 4.  So  we  must  have  that  F is  convex  on  X which  means 

n 

also  that  X must  be  convex.  Let  D(0)  e A and  P be  a permutation  matrix  ( 6^_.+^) ", 


(n+1  = 1) . Then 


(2.26) 


1 n i n 

r-  l PTXBHP  ) = ( l 0.  /n)  I . 

i=l  i=i  1 


T — 

Here  by  P we  denote  the  transpose  of  P.  Therefore  if  0 e X then  0 = (b, . . . ,b)  £ X 


(b  » 1 0.  /n)  . This  in  particular  implies  that 

i-1 


(2.27)  if  0 e x,  a < 0,  then  a £ X. 

P^inition  2.1.  Let  X c.  F>  . The  set  X is  called  strongly  convex  if  X is  convex  and  the 
condition  (2.27)  is  satisfied. 

Theorem  2.2.  Let  X be  a strongly  convex  set  in  f"  which  contains  at  least  one  point  a , 

(2.28)  a,  > a.  >...>  a 

12  n 

Let  F : X -*•  F . Assume  that  F £ (X)  . Consider  a spectral  function  if  : A + F (A  c Hn> 

where  if  and  A ;re  given  by  (2.11)  and  (2.25)  accordingly . Then  if  is  convex  on  A if  and 


only  if  F is  convex  on  X and 


3x7  <a>  - 3x7  <a>  i ta-  <a) 

12  n 


for  any  a £ X.  Moreover , 4>  is  strictly  convex  on  A,  i.e. 


(2.30)  if  (cA  + (l-c)B)  < cif(A)  + (l-c)if(B),  A*B,  0<c<l  , 

if  and  only  if  F is  strictly  convex  on  X and 

(2. 3D  ^7  (a)  > -3£-  („)  if  « > a . 

i l+l 

To  prove  the  theorem  we  need  the  following  theorem  of  Ostrowski  [11]  (Theorems  VII  and 


I 


f 


Theorem  2.3.  Let  X and  F satisfy  the  assumptions  of  Theorem  2.2.  Then  F satisfies  (2.29) 
if  and  only  if 

(2.32)  F(a)  < F(B)  if  a < 6. 

Moreover 

(2.33)  F(a)  < F(B)  if  a < B and  a B 
if  and  only  if  the  condition  (2.31)  holds . 

Proof.  Assume  first  that  F is  convex  on  X.  So  if  1(A),  A(b)  e X then 

(2.34)  F(— j (F(A  (A) ) + F(A(B)))  . 


According  to  Theorem  2.3,  the  assumption  (2.29)  implies 


(2.35) 


F(M^BL)  < F( 


A (A)  -f  A (B) . 
2 ' 


by  the  virtue  of  (2.15).  This  shows  that  $ is  convex  on  A.  Assume  furthermore  that  F is 
strictly  convex  on  X.  So  if  A (A)  / A(B)  the  inequality  sign  holds  in  (2.34).  This  implies 
(2.30).  Suppose  that  A (A)  = A(B)  but  A ^ B.  According  to  Lenma  2.1  A (A+B)  ^ (A(A)  + A(B)) 
So  the  additional  assumption  (2.31)  yields  the  inequality  sign  in  (2.35)  according  to  Theorem 
2.3.  This  manifests  that  <p  is  strictly  convex  on  A.  Assume  now  that  <(>  is  convex  on  A. 
In  particular  $ is  convex  on  n A.  This  immediately  implies  that  F is  convex  on  X. 
Furthermore  if  $ is  strictly  convex  then  F is  strictly  convex.  Let  B e X.  So  D(8)  e A. 
Assume  that  a < B.  Then  D(a)  e A.  The  classical  result  of  [7, sec.  2.19]  states  that 

(2.36)  G6  = a , 

where  G is  some  doubly  stochastic  matrix.  The  Birkhoff  theorem  implies 

k k 

(2.37)  G = l a P , a > 0,  l a = 1 

i=l  i=l 


and  p^  is  a permutation  matrix.  So 

(2.38)  D(a) 


l aiP.D(B)P^  . 


So  the  convexity  of  0 implies 
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3 . some  characterization  of  the  spectral  radius 

Let  A be  an  n x n non-negative  matrix  such  that  there  exists  two  positive  vectors  u,  v 
satisfying 

(3.1)  Au  = r(A)u,  AvT=  r (A) V , uT  = (u  ,...,u  ) > 0,  vT  = (v  , . . . ,v  ) > 0 

In  In 


Assume  the  normalization 
(3.2) 


7 u.v.  = 1 . 

i=l 


Let  be  the  set  of  probability  vectors 


(3.3) 


Pn  - {al°  = (“l an} ' “i  i°'  J,  “i  = X)  • 

i=l 


In  (5,  Sec.  3]  it  was  manifested 

Theorem  3.1.  Let  A be  an  n * n non-negative  irreducible  matrix  having  positive  entries  on 

the  diagonal  (or  fully  indecomposable , see  Remark  3.3  ill  [5]).  Then  for  any  a e P , with 

n (Ax) . 

positive  entries  (cu  > 0).  The  function  f(x)  = 7 log  — “ has  a unique  critical  point 

i=l  i 

E = (£.,...,£  ) in  the  interior  point  of  P (£.  > 0)  which  must  satisfy 
1 n — n l 

(Ax)  j n (A£)  . 


(3.4) 


min  T a.  log = 7 a.  log  — . 

. i x.  . L . l £. 

x>0  i=l  i i=l  l 


Thus,  if  a is  chosen  to  be 
(3.5) 


a = (u.v  ,...,u  v ) , 

11  n n 


where  u and  v satisfy  (3.1)  - (3.2)  then 

n 

l 

i=l 


(3.6) 


(Ax)  i 

l u.v.  log  — log  r(A) 


since  x = u is  a critical  point  of  f(x) . 

Prom  Theorem  3.1  we  get 

Theorem  3.2.  Let  A be  an  n » n non-negative  matrix  such  that  r(A)  >0.  Then 

n (Ax) ^ 

(3.7) 


sup  inf  7 ^ ~ = log  r (A) 


aeP  x>0  i=l 
n 
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1 1 


Suppose  that  there  exists  a positive  vector  u satisfying  (3.1).  Assume  that 


(3.8) 


n (Ax) . 

inf  1 a.  log - 

x>0  i=l  1 xi 


log  r(A) 


Then  the  vector  v 


(3.9) 


v = (“i^i w 


I 


I ■ 


' 

1 

1 

I 

i 

» 1 


I 


! 

■ ( 

E:| 


fulfills  (3.1).  In  particular  if  A is  irreducible  then  a is  unique  and  given  by  (3.5). 

Proof : As  the  left-hand  side  of  (3.7)  is  a continuous  function  of  A it  is  enough  to  prove 

(3.7)  for  A positive.  Let  u > 0 be  the  corresponding  eigenvector  of  A.  So 

n (Ax) ^ n (Au) ^ 

inf  l a.  log  — - 1 £ l a.  log  = log  r(A) 

x>0  i=l  1 Xi  i=l  1 Ui 


for  any  a such  that  ][  = 1. 


Thus 


i=l 


(Ax), 


inf  £ a . log  — log  r (A) 

...n  I ^ J 


sup 

a£P  x>0  i=l 
n 


The  above  inequality  together  with  (3.6)  yields  (3.7).  Suppose  that  (3.8)  holds.  If  u > 0 


n (Ax) ^ 


satisfies  (3.1)  then  x = u is  a minimal  point  for  f(x)  = £ a.  log 


So 


i=l 


3f 

3x . 
3 


n a. a.  . a. 

= V -JLil  _ _1 

(Ax).  x. 

X=U  1=1  1 1 


-1 


= r(A)  1 V a.u.^a, . - a.u. 
;=u  i=l  1 1 j 


This  shows  that  v given  by  (3.9)  is  a left  eigenvector  of  A corresponding  to  r(A) . If  A 
is  irreducible,  then  u and  v are  unique  up  to  a multiple  of  a positive  scalar.  Thus  a is 
of  the  form  (3.5)  and  since  a £ P^ , a unique.  The  proof  of  the  theorem  is  completed. 

We  now  bring  an  extended  version  of  Theorem  3.2  which  includes  (3.7)  and  the  Donsker- 
Varadhan  characterization  [2]  as  its  special  cases. 

Theorem  3.3.  Let  ¥ : F ->•  3R  be  a continuous  convex  function  on  F . Define  $ : F+  -»  IR 

(3.10)  <Mx)  = ¥(log  x)  . 

Let  A be  an  n x n non-negative  matrix  such  that  r(A)  > 0.  Assume 

(3.11)  ¥'  (log  r (A) ) > 0 . 

-10- 





. . 


sup  inf 

cuP  *>0 
n 


n /(Ax)  A 

i °i*br) 


♦ (r (A) ) . 


Assume  that  the  inequality  sign  holds  in  (3.11)  and  suppose  that  there  exists  a positive  vector 


(3.1).  If 


n /(Ax)  A 


♦<r(A))  , 


then  the  vector  v (3.9)  satisfies  (3.1).  In  particular  if  A is  irreducible  then  a is 


unique  and  given  by  (3.5). 


Proof.  Let  t.  = log  r(A) , ¥• (t _)  - e.  Then  the  convexity  of  ¥ implies 


y(t)  > ¥(tQ)  + (t-t0>r  <t0)  . 


n / (Ax)  A 

l ”‘*br)  * 


(*(r (A) ) - e log  r (A) 


n (Ax)  i 

+ e l a log  — , a € P 

i-1  1 i 


As  e >_  0 from  Theorem  3.2  and  the  above  inequality  we  get 


sup  inf 

aeP  x>0 
n 


n /(Ax).\ 

& *■ 


*(r(A) ) . 


Since  t is  continuous  we  may  assume  that  A is  positive.  By  choosing  x - u the  left-hand 
side  of  (3.15)  we  deduce  an  opposite  inequality  of  (3.15) . This  establishes  (3.12) . In  case 
the  e > 0 we  use  the  arguments  of  Theorem  3.2  to  analyze  the  equality  (3.12).  End  of  proof. 

Letting  ¥<x)  “ ex  in  Theorem  3.3  we  obtain  the  Donsker-Varadhan  characterization  [2]. 
Corollary  3.1.  Let  the  assvsiiptions  of  Theorem  3.2  hold.  Then 


(3.16) 


n (Ax) ^ 

sup  inf  £ a. - r(A) 

a«P  x>0  i-1  Xi 

n 


Suppose  that 


-11- 


inf  sup  I a 7 - -777  . 

oeP  x>0  i=l  1 (Ax,i  r(A> 
n 


Assume  that  there  exists  a positive  vector  u satisfying  (3.1)  and  su 
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and  B is  reducible  if  and  only  if  A is  reducible  (e.g.  [8,  chap.  8].  Again,  as  in  the 
proof  of  Theorem  3.2  one  may  assume  that  B is  positive.  By  letting  x to  be  equal  to  the 
positive  eigenvector  u of  A we  immediately  deduce 


(3.23) 


l «i 


in„f  .n  ,L,  “i  CAX)  . - r(A) 
a«P  x>0  i-l  i 

n 


Let  a be  given  by  (3.5).  Obviously  for  any  x > 0 and  y - Ax 

n x;  n (A  ^y) , n (By)  . 


(3.24) 


it  a.  it  \n  ^ »» 

ja  Vi  - j,  Vi  -Pf-  - ' - ,1  Vi  — 


From  Corollary  3.1  it  follows 


(3.25) 

So 

(3.26) 


(By)i 

l Vi  — 1 r(B) 


i-1 


? Vi  735o"- r • r(B)  = 7Ti 


i-l 


(A) 


and  the  equality  sign  holds  if  x = u.  This  establishes  (3.20) . The  equality  (3.21)  is 
analyzed  in  the  same  way  as  in  Theorem  3.2. 

Remark  3.1.  Theorem  3.4  does  not  hold  for  arbitrary  non-negative  matrices,  take  for  example 
A to  be  a permutation  matrix  P ft  I.  Therefore  Theorem  3.4  is  not  a special  case  of  Theorem 
3.3. 
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4.  Convexity  properties  of  the  spectral  radius 


Let  A be  an  n * n non-negative  matrix.  Consider  the  matrix  A + D,  D e D . Assume 

n 


that  the  eigenvalues  of  A D arranged  in  the  order 


Re  A (A)  > Re  A, (A)  >...>  Re  A (A) 
i z — — n 


p(D)  = Ax(A+D) 


We  claim  that  p(D)  is  real.  If  d is  non-negative  this  fact  is  a consequence  of  the  Perron 

Prober ius  theorem.  For  an  arbitrary  D consider  A + D + al 

<4-3>  \ (A+D+al)  = Ak(A+D)  + a,  k = l,...,n 

for  a big  enough  A + D + al  _>  0 and  (4.3)  implies  that  p(D)  is  real.  Moreover  by  consider 

ing  the  matrix  B = A + D + al  and  using  the  Donsker-Varadhan  characterization  for  B we  get 

the  following  characterization  for  p(D) 

<4-4)  p(D)  = sud  L (D,a) 

“£Pn 


Here  L..(D,a)  is  a linear  functional  on  V 
1 n 


n (Ax) . 


VD'“>  = l “idi  + l “i  ~T~  ' 

1=1  X>0  1=1  i 


o = (a. , . . . ,a  ) , D = diag{d  , . . . ,d  } . 

n in 

It  is  a standard  fact  (4.4)  and  (4.5)  inqply  the  convexity  of  p(D)  on  the  set  D . More  pre- 

n 

cisely  we  have: 

Theorem  4.1.  Let  A be  a fixed  n x n non-negative  matrix.  Assume  that  p(D) . D e D , is 


iven  by  (4.2) . Then  p(D)  is  a real  valued  convex  functional  on  D 


l4-6>  pUDj+Dj)^)  _<  (plD^  + p(D2))/2  . 

Moreover  if  A is  irreducible  then  the  equality  sign  holds  in  (4.6)  if  and  only  if 


I I 

t Ir 


Proof . As  we  pointed  out  (4.6)  is  a consequence  of  (4.4).  So  it  is  enough  to  analyze  the 
equality  case.  Let 

(4.8)  Ax  - A + (Dj+Dj)/^,  AjU  • r^U,  A^V  - ^v,  ^ /2)  . 

As  A is  irreducible  we  may  assume  that  u,  v > 0 and  the  normalization  (3.2)  holds.  Let  a 
be  given  by  (3.5).  So 

n (Ax)  n (A  u) 

(4.9)  L <(D  +D  )/2,a)  = inf  l a = l a — ± - r.  . 

x>0  i«l  1 Xi  i-1  1 Ui  1 


b 


If  we  apply  the  results  of  Section  3 in  [5J 

n (Bx) 

(4.10)  f(x,B)  = l a — -i 

i-1  1 xi 

where  B + bl  is  irreducible  matrix,  for  some  positive  b,  then  f(x,B)  has  a unique  critical 
point  in  the  interior  of  P^  which  must  be  the  minimum  point  (f(x)  - +®  on  the  boundary  of 
P ) . The  equality  sign  in  (4.6)  implies 

(4.11)  L1(D1,a)  = p (D1>  , L1(D2,a)  - p(D2)  . 


(4.12)  ftx^+D^  ^ ffu.A+D.^)  = p (Dj.)  , f(x,A+D2>  >_f(u,A+D2>  - p(D2> 

The  uniqueness  of  the  minimal  point  of  f(X,B)  implies 

(4.13)  (A+Dj^u  = ptDj^U,  (A+D2)u  - p(D2)u  . 

As  u > 0 (4.7)  follows  the  above  equality.  The  proof  of  the  theorem  is  completed. 

The  inequality  (4.6)  extends  Cohen's  result  (1).  Let  A be  a non-negative  matrix  such 

that  r (A)  > 0.  Clearly,  for  any  D e D , r(e°A)  is  also  positive.  Define 

n 

(4.14)  R(D)  - log  r(e°A) 


According  to  Theorem  3.2, 


R(D)  - sup  L (D,a)  , 

aeP 

n 


" MS  /• 
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(4.16)  L (D,a)  = l a d.  + inf  £ a.  log  . 

i-1  1 1 x>0  i-1  1 xi 

Combining  (4.15)  and  (4.16)  and  using  the  uniqueness  result  stated  in  Theorem  3.1  as  in  the 
proof  of  Theorem  4.1  we  deduce. 

Theorem  4.2.  Let  A be  a fixed  n x n non-negative  matrix  having  a positive  spectral  radius 


Assume  that  R(D)  is  given  bv  (4.14).  Then  R(D)  is  a convex  functional  on  D 


(4.17)  R((D1+D2)/2)  £ (RfD^  +R(D2))/2  . 

Moreover  if  A is  irreducible  and  the  diagonal  entries  of  A are  positive  (or  A is  fully  in 


decomposable)  then  the  eauality  sign  holds  in  (4.17)  if  and  only  if  (4.7)  holds  for  some  a 


Assume  that  A,  B e H and  furthermore  A is  positive  definite  ((Ax,x)  > 0 for  x ^ 0) 
n 

1/2  1/2 

BA  is  similar  to  A BA  . This  shows  that  X, (BA)  is  a convex  functional  on  H for 

1 n 

a fixed  positive  definite  A.  If  in  addition  A has  non-negative  entries  then  X^DA)  is  con 

vex  on  This  result  does  not  apply  in  general  for  non-negative  matrices.  For  example,  take 

A to  be  a permutation  matrix  P f I.  However,  X, (DA)  is  convex  on  D*  - the  set  of  n » n 

i-  n 

non-negative  diagonal  matrices  if  A-1  is  an  M-matrix. 

Theorem  4.3.  Let  A 1 be  an  M-matrix.  Then  r (DA)  is  a convex  functional  on  d\ 


/(d1+D2>  \ 1 

2 Aj  < j (r  (DjA)  + r (D2A) ) . 


Moreover  if  A is  irreducible  then  the  equality  sign  in  (4.18)  holds  if  and  only  if 


sitive  a 


Proof.  Using  the  continuity  argument  we  may  assume  that  in  the  decomposition  (3.22)  B is 
positive  (irreducible),  i.e.  A is  positive  (irreducible).  Thus  if  all  diagonal  elements  of 
Dq  » diag  {d°,...,d°}  are  positive  then  DQA  is  positive  (irreducible).  According  to  the 
Perron-Froberius  theorem  r (DQA)  is  a simple  root  of  the  det(XI  - DQA)  = 0.  By  the  implicit 
function  theorem  r (DA)  is  an  analytic  function  of  D in  the  neighborhood  of  DQ.  Then  the 
convexity  of  r (DA)  would  follow  if  we  show  that 


I 


(4.20) 


r(DA)  > r(D0A)  + l (d^J0*)  Q , 

i-1 


for  any  DQ  with  positive  diagonal  elements.  Let  £.1  be  the  eigenvectors  corresponding  to 


DqA  and  A Dq 


(4.21) 


DqA£  - r(DQA)C,  A DQn  - r (DQA)  tl  , 


o*£=  0 < n « (tijf  •••'”„> * l * 1 


It  can  be  shown  that 


(4.22) 


ar(DA) 


3d . 


T 3n  ^i^i 

0 " n 3d-  “ r(D0A)  .(0)  i “ 1"*"n 

D 1 d^ 


This  can  be  done  by  bringing  DQA  to  the  Jordan  form  and  using  the  simplicity  of  r(DQA) . Sei 


for  example  [10,  II,  55.4).  Thus  (4.20)  is  equivalent  to 

n d. 


(4.23) 


r (DA)  ir(D0A)  I -^5^ 
i-1  d± 


This  inequality  was  established  in  15).  It  follows  directly  from  (3.26).  Indeed  suppose 
that  D has  positive  diagonal  elements  and  let 


(4.24) 

Then  according  to  (3.26) 


DAw  = r(DA)w,  w = (w,  ,...,w  ) > 0 
l n 


JL  i 

r(D0A)  — i i (DgAw)^ 

n *^i  wi  1 n “i 

X 5ini  7(0 T (DAw)  . = r (DA)  eini  .(0)  ' 

i=l  a.  i i=l  d. 


which  establishes  (4.23)  for  D with  positive  diagonal.  So  (4.18)  holds  in  the  interior  of 


D+.  The  continuity  argument  implies  the  validity  of  (4.18)  on  D* . Suppose  that  A is  also 
n n 


irreducible.  Then  B in  the  decomposition  (3.22)  is  also  irreducible,  since  the  inverse  of 
block  triangular  matrix  is  also  a block  triangular  one.  As  in  the  proof  of  Theorem  4.1  strict 
inequality  holds  in  (4.23)  unless  D„A  and  DA  have  the  same  positive  eigenvector.  So 


D - aDQ  for  sane  a > 0.  This  shows  that  we  have  strict  inequality  in  (4.18)  unless  (4.19) 


holds  provided  that  DQ  (which  is  either  or  D2)  have  positive  diagonal.  The  proof  of 


the  theorem  is  completed. 


We  conclude  this  section  by  pointing  out  that  the  convexity  of  r (DA)  on  D+  is  a strong- 


er result  then  the  convextiy  of  log  r(eCA)  on  D . indeed,  let 

n 


(4.25) 


D0  - e °,  CQ  - {q<°> qj°>},  q<°>  = log  d{0),  i - 1 . 


Suppose  that  log  r(e®A)  is  convex  at  Q = QQ.  This  means 


(4.26)  log  r(eCA)  > log  r(DQA)  + r^A)"1  | (q^0* ) , Q - diagt^ qj 


As  in  the  proof  of  Theorem  4.3 


(4.27) 


9r (e^A)  T 3e°  . , , 

3q,  Q~Qn  " n 9q  Qn  5 “ r(D0A)ni5i'  i ~ 1"”'n 


where  n,C  given  by  (4.21), 


Thus  (4.26)  is  equivalent 


(4.28)  r (e 


2 4qA« ini  n / di  \ iT1i 

3A)  > r(D0A)n  — j =r(D0A)  n l-±\  ,q. 

W 1“1\d1  > 


log  d,  , i — 1, . . . ,n  . 


Using  the  relation  between  the  arithmetic  and  the  geometric  means  from  (4.23)  we  get 


(4.29) 


n qi 


I'1"1 


rlA’ 1 jfer  «i»i  i *'v>  Afc 


That  is  the  convexity  of  r (DA)  at  dq  c C*  implies  the  convexity  of  log  r(eSA)  at 


PQ  ” lo9  Dq.  This  demonstrates  that  the  convexity  of  r (DA)  on  D*  implies  the  convexity  of 
log  r(e  A)  on  D . On  the  other  hand  if  A is  a permutation  matrix  ^ I then  r(DP)  is  not 


convex  on  D (for  details  see  [5],  Section  3). 


r 


5.  Convex  functions  on  the  spectrum  of  totally  positive  matrices 

A real  valued  n * n matrix  is  called  a totally  (strictly  totally)  positive  matrix  of 
order  k if  all  minors  of  A of  order  less  or  equal  to  k are  non-negative  (positive) . We 
denote  these  matrices  by  TP^(STP^).  For  j = n we  call  these  matrices  simply  by  TP(STP) . A 
matrix  A is  called  oscillating  if  A is  TP  and  some  power  of  A is  STP.  It  is  known  that 
a TP  matrix  if  oscillating  if  and  only  if 


(5.1) 


aii  ” °*  ai(i-l)  ” °'  ai(i+l)  > 0.  i-l. 


,n,  A * L 0 • 


In  that  case  A is  totally  indecomposable. 


If  A is  TPj  then 


(5.2) 


X1(A)  > X2(A)  >...>  Aj  (A)  > |Xfc(A)  |,  k - j + If • • • #n  • 


If  A is  STPj  then  we  have  strict  inequalities  in  (5.2).  See  [6)  and  [9)  for  proofs  of  these 

results  and  more  properties  of  these  matrices.  Let  A be  TP..  Define  $ : A -+  it  (^  c D ) 

3 n 

as  follows 

(5.3)  4(D)  - F(log  X (e°A) , . . . , log  X.(e°A))  . 


As  in  Section  2 w"  were  looking  for  necessary  and  sufficient  conditions  on  F which  imply 

that  ♦ is  a convex  function  on  A c D for  any  A which  is  TP_. . It  turns  out  that  we  have 

an  analogous  result  to  Theorem  2.2.  To  do  so  we  need  few  notations  and  definitions.  Let 

a « (Hj,,..,o.)  and  6 ■ (6^,...,Bj)  and  j < n.  We  define  a « B if  (2.12)  holds  for 

k « l,...,j.  Thus  if  a « (a,,..., a ),  8 « (6,,...,B  ) and  a < B then  a <<  B.  Conversely, 

In  In 

if  a « B we  can  extend  a to  o and  B to  8 such  that  a < B.  A set  X c is 

called  a super  convex  if  X is  convex  and 


(5.4) 


if  B £ X,  a « B,  then  o t X 


Clearly  X is  super  convex  in  if  and  only  if  it  could  be  extended  to  X £ r"  such  that 

X is  strongly  convex  in  It"  . Using  the  above  arguments  and  Ostrowski's  result  (Theorem  2.3) 
we  get 

Lemma  5.1.  Let  X be  a super  convex  set  in  . let  F s X -*■  It  . Assume  that  F e C*1'  (X)  . 
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Then 


if  and  only  if 


F(oO  < F(B)  if  a « B 


for  any  a e X.  Moreover  strict  inequality  in  (5.5)  holds  for  a ft  6 if  and  only  if 

(5-7)  ir  <a>  3^7  (“>  — ai > vr  ^7  (“>  * 0 if  aj  * 0 • 

Assume  that  A is  TP^.  Denote 

(5.8)  X(j)(A)  - (X1(A),...fXj(A)),  log  X(j)(A)  - (log  X^A)  , . . . ,log  X..(A))  . 

Theorem  5.1.  Let  A be  an  n x n non-singular  TP^  matrix.  If  i < n then 

(5.9)  log  X(5,(e  1 A)  « ± log  x'3’  (e  XA)  + log  X(j)  (e  2A) ) . 


If  j « n then 


(5.10) 


(d1+d2)/2  dx  d2 

log  X(e  A)  < — [log  X(e  A)  + log  X(e  A)) 


If  in  addition  A satisfies  (5.1),  or  more  generally  A is  totally  indecomposable,  then 


(5.11) 


(1)  1 M\  D1  (4\  ^2 

log  X°'  (e  A)  - | [log  X(3)  (e  XA)  + log  X(j)  (e  2A)  ] 


for  any  1 < j < n if  and  only  if  (4.7)  is  satisfied  for  some  a. 
Proof.  Denote  by  C^(A)  the  k-th  compound  of  A.  Thus 


(5.12) 


D VD) 

Ck(eu)  - e K 


where  * is  well  defined  map  * : D -+D  . It  is  easy  to  see  using  the  properties  of  the 

* k n .n. 

'fc'  rst 

compound  matrices  that  ^ is  a linear  map.  According  to  Theorem  4.2  log  r(e  C^fA) ) is  con- 
vex on  D for  k = l,...,j.  Note  that  the  non-singularity  of  A implies  that  r(C  (A) ) > 0. 
<k> 

\ (D) 

Thus  log  r(e  C.(A) ) is  convex  on  V . Let 
k n 


isable  then 


is  strictly  convex  if 


and  only  if  F is  strictly  convex  and  satisfies  (5.7).  In  case  j - n,  4>  is  convex  (strictly 
convex  provided  that  A is  totally  indecomposable)  if  and  only  if  F satisfies  the  assumptions 


mk 


It  is  well  known  that 


\ (D)  - l log  X.  (e  a)  . 
K i-1  1 


^(D)  - log  r(Ck(«DA>) 


Therefore  (D)  is  convex  on  for  k - l,...,j.  This  is  equivalent  to  (5.9)  for  j < n. 


(5.15)  R (D)  - log  det(eDA)  » £ d.  ♦ log  det(A) 

i-1  1 

This  verifies  (5.10)  if  A is  a TP  matrix.  Suppose  that  in  addition  A is  totally  indecom- 
posable. According  to  Theorem  4.2  we  have  a strict  inequality  in  (4.17)  unless  (4.7)  holds. 
Thus  (5.11)  can  be  satisfied  if  only  (4.7)  holds.  Trivially  (4.7)  implies  (5.11).  The  proof 
of  the  theorem  is  completed. 

Theorem  5.2.  Let  X be  a super  convex  set  in  for  1 ^ j < n (a  strongly  convex  con- 
taining a point  a,  >...>  an>  if  j = n)  . Let  P i X +*.  Assume  that  F c C(1>  (X)  . Let 

A be  a given  n x n non-singular  TP..  matrix.  Consider  a spectral  function  ♦ : A * » , 
given  by  (5.3) , where  A is  a convex  set  in  D such  that 

(5.16)  log  X^5(eDA)  c x,  D € A . 


Then,  for  all  such  A,  $ is  convex  if  and  only  if  F is  convex  on  X and  satisfies  (5.6)  in 


of  Theorem  2.2. 

Proof.  A proof  of  this  theorem  can  be  achieved  by  modifying  in  the  obvious  way  the  proof  of 
Theorem  2.2.  In  fact,  all  the  arguments  of  the  proof  of  Theorem  2.2  carry  over  if  one  notices 
that  the  identity  matrix  is  TP. 
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20  ABST'IACT  (Continue  on  reveree  aide  If  neceaamry  end  Identity  by  block  number) 

In  this  paper  we  characterize  all  convex  functionals  defined  on  certain 
convex  sets  of  hermitian  matrices  and  which  depend  only  on  the  eigenvalues  of 
matrices.  We  extend  these  results  to  certain  classes  of  non-negative  matrices. 
This  is  done  by  formulating  some  new  characteriz-ations  for  the  spectral  radius 
of  non-negative  matrices,  which  are  of  independent  interest. 
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